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ABSTRACT 


The  goal  of  this  thesis  is  to  investigate  the  initial  search  planning  phase  of  the  Search 
and  Localization  Tactical  Decision  Aid  (SALT)  developed  by  MET  RON,  Incorporated 
of  McLean,  Virginia.  SALT  is  a  Computer  Assisted  Search  (CAS)  program  intended  for 
use  by  P3  UPDATE  IV  crews  to  assist  them  in  optimal  deployment  of  a  sonobuoy  field 
to  prosecute  a  submarine  threat. 

The  initial  search  planning  phase  of  SALT  takes  as  user  inputs  environmental  data, 
an  initial  elliptical  Search  Probability  Area  (SPA),  an  assumed  target  motion  model,  and 
the  duration  of  the  search.  Outputs  include  a  recommended  sonobuoy  pattern  and  the 
probability  of  detection  of  this  pattern.  The  investigation  of  this  phase  of  the  algorithm 
is  conducted  in  two  parts,  first,  a  series  of  simulation  routines  is  used  to  ensure  that  the 
probability  of  detection  of  the  sonobuoy  patterns  generated  by  SALT  is  mathematically 
correct.  Second,  these  same  simulation  routines  are  used  to  determine  if  there  arc  al¬ 
ternate  sonobuoy  patterns  that  result  in  higher  probabilities  of  detection. 
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I.  INTRODUCTION 


Tactical  Decision  Aids  (TDA's)  have  been  developed  and  implemented  on  comput¬ 
ers  in  various  forms  to  assist  users  in  making  decisions  on  a  wide  variety  of  topics.  One 
prominent  area  where  TDA's  have  proven  to  be  of  value  is  in  the  area  of  Antisubmarine 
Warfare  (ASW).  In  the  search  theory  phase  of  the  ASW  problem,  in  particular,  there 
is  an  extremely  large  amount  of  data  which  must  be  assimilated  and  presented  to  the 
user  in  a  manner  conducive  to  rapid  evaluation  so  that  a  sound  decision  can  be  made 
and  acted  upon.  Several  TDA's  have  been  developed  to  aid  the  user  in  the  area  of  search 
theory.  These  TDA's  have  been  grouped  together  under  the  title  Computer  Assisted 
Search  (CAS).  A  common  thread  which  runs  through  all  of  the  CAS  programs  is  the 
concept  of  a  "probability  map".  It  is  this  probability  map  which  allows  all  of  the  input 
data,  mathematical  calculations,  and  time  step  updates  to  be  output  as  a  2-dimcnsional 
graphic  display  that  is  easily  understood  by  the  user  in  terms  of  relative  probabilities  on 
the  location  of  the  target  of  interest. 

Many  CAS  programs  represent  the  stochastic  aspect  of  target  motion  by  a  Monte 
Carlo  simulation  of  a  large  number  of  possible  target  tracks.  As  information,  cither 
positive  (detection)  or  negative  (no  detection  after  a  period  of  time),  is  obtained  from  the 
seal ch,  the  probabilities  of  the  individual  tracks  are  updated  and  displaced  as  a  proba¬ 
bility  map  at  the  user  requested  time.  The  most  prominent  TDA's  of  this  type  arc 
VPCAS  and  PACSIiARCI  I.  This  type  of  TDA  possesses  an  inherent  simplicity  which 
piomotcs  ease  of  understanding  of  the  underlying  concepts.  The  use  of  the  Theorem  of 
'I otal  Probability  and  Bayes'  Thcoiem  arc  easily  applied  to  small  scale  problems  which 
can  be  worked  by  hand  and  the  results  naturally  generalize  to  large  scale  Computer  As¬ 
sisted  Seal  cli  problems.  The  difficulty  in  this  type  of  approach  is  that  the  specific  tracks 
selected  for  the  probability  map  do  not  cover  all  of  the  possibilities  of  target  motion. 
[Ref.  1:  pp.  1I-1.1I-2J 

In  an  attempt  to  negate  the  above  difficulty  with  the  Monte  Cailo  approach  to 
target  motion,  a  new  algorithm  entitled  Search  and  Localization  Tactical  Decision  Aid 
(SALT)  has  been  developed  by  ML!  RON,  Incorporated  of  McLean.  Virginia.  SALT 
is  a  CAS  algorithm  intended  for  use  by  P3  L’PDATL  IV  crews  to  assist  them  in  the  op¬ 
timal  deployment  of  a  sonobuoy  field  to  prosecute  a  submarine  threat.  The  SALT  al- 
goiithm  uses  a  discrete  time  and  space  Markov  Chain  to  simulate  target  motion.  The 
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state  space  of  this  Markov  Chain  is  a  five  element  vector  (position,  velocity,  depth).  The 
number  of  discrete  cells  in  the  state  space  is  21  entries  each  for  latitude  and  longitude, 
eight  entries  each  for  north-south  and  cast-west  velocity  components,  and  two  entries 
for  depth.  Thus  the  current  state  of  the  target  is  represented  probabilistically  as  56,44S 
five-dimensional  cells.  This  analuic  approach  to  solving  the  problem  of  unknown  target 
motion  allows  for  many  more  possibilities  of  target  motion  than  the  simplistic  track 
technique  used  in  previous  CAS  models.  [Ref.  1:  pp.  11-46,  11-47] 

The  initial  search  planning  phase  of  SALT  takes  as  user  inputs  environmental  da»a, 
an  elliptical  Search  Probability  Area  (SPA),  an  assumed  target  motion  model,  and  du¬ 
ration  of  search.  The  algorithm  then  outputs  a  recommended  sonobuoy  deployment 
pattern  along  with  the  probability  of  detection  of  this  pattern.  [Ref.  2:  p.  1]  This  initial 
search  planning  phase  is  the  topic  which  will  be  investigated  here.  The  investigation 
consists  of  two  parts,  first,  the  sonobuoy  patterns  recommended  by  SALT  v  ill  be  tested 
using  FORTRAN  simulations  to  ensure  that  the  probability  of  detection  for  each  pat¬ 
tern  as  output  by  SALT  is  mathematically  accurate.  Second,  alternative  sonobuoy  pat¬ 
terns  will  be  explored  with  the  goal  of  developing  patterns  that  improve  on  the  SALT 
probability  of  detection. 

We  shall  begin  in  Chapter  11  by  discussing  the  three  types  of  motion  models  used 
by  SALT  to  develop  initial  probability  maps.  Chapter  111  follows  with  a  presentation 
of  how  buoy  patterns  arc  developed  and  evaluated  based  on  the  input  data  for  each  type 
of  target  motion.  Chapter  IV  provides  the  documentation  of  the  FORTRAN  simu¬ 
lation1:  that  are  used  to  verify  the  probability  of  detection  for  the  SALT  generated  pat¬ 
terns.  In  Chapter  V.  the  results  obtained  from  the  FORI  RAN  simulation  will  be 
presented  and  compared  to  the  output  from  the  SALT  algorithm.  Chapter  VI  explores 
the  use  of  alternative  sonobuoy  patterns  with  the  goal  of  developing  put'erns  which  will 
improve  the  probability  of  detection  of  those  patterns  currently  used  in  the  SALT  pro¬ 
gram.  Chapter  VII  is  a  summary  of  the  results.  The  inputs  to  the  SALT  algorithm  used 
to  obtain  the  sonobuoy  patterns  for  testing  are  listed  in  Appendix  A.  'I  he  source  code 
for  examples  of  the  FORTRAN  simulation  routines  is  presented  in  Appendix  B. 


II.  SALT  MOTION  MODELS 

There  are  three  different  models  of  target  motion  built  in  to  the  SALT  algorithm. 
These  models  are  a  patrolling  target,  a  transiting  target,  and  a  fleeing  datum  target.  The 

specifics  of  each  type  of  target  motion  are  discussed  below.  In  general,  each  scenario 

\ 

begins  with  an  elliptical  SPA  which  is  input  by  the  user.  To  specify  a  SPA  requires  seven 
input  parameters:  date-time-group,  latitude  and  longitude  of  the  SPA  center,  semi-major 
axis  (nm),  semi-minor  axis  (nm),  orientation  (deg  N),  and  containment  percent.  This 
SPA  represents  a  Bivariate  Normal  distribution  for  the  target's  location  that  has  been 
obtained  from  an  external  source.  [Ref.  3:  pp.  12,  21] 

A.  PATROLLING  TARGET  MODEL 

In  the  patrolling  target  model,  the  target  is  assumed  to  be  randomly  patrolling  in  the 
vicinity  of  the  SPA.  The  target  selects  an  initial  course  from  a  Uniform  (0,  360)  dis¬ 
tribution  and  an  initial  speed  from  a  truncated  triangular  distribution.  This  triangular 
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distribution  takes  as  input  the  average  target  speed.  The  distribution  ranges  from  half 
of  the  input  speed  to  two  times  the  input  speed  with  the  average  speed  being  eight  times 
as  likely  as  cither  of  the  endpoints.  These  distributions  are  shown  in  Figure  1  below. 


figure  1.  Patrolling  Target  Course  and  Speed  Distributions 
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The  target  is  assumed  to  continue  on  this  course  and  speed  for  an  amount  of  time 
determined  from  an  exponential  distribution  with  mean  time  between  course  changes 
being  input  by  the  user.  At  the  completion  of  this  leg  in  the  target's  track,  an  inde¬ 
pendent  course,  speed,  and  time  on  leg  are  selected  from  the  same  distributions  and  the 
process  continues  throughout  the  scenario.  [Ref.  3:  pp.  10-11,20] 

The  values  of  the  input  parameters  for  a  patrolling  target  were  set  as  follows  to  ob¬ 
tain  SALT  output  for  testing: 

•  Average  target  speed  (knots)  4 

•  Average  time  between  course  changes  (hours)  3 

B.  TRANSITING  TARGET  MODEL 

The  transiting  target  motion  model  assumes  that  the  target  has  a  preferred  course 
and  speed  of  advance  (SOA)  to  travel.  In  this  model,  the  user  is  required  to  input  five 
parameters.  The  parameters  are  average  target  course  (deg  N),  maximum  variation  in 
target  course,  average  SOA  (kts),  maximum  variation  in  SOA,  and  average  time  between 
course  changes.  For  each  leg.  the  target  selects  a  course  from  a  truncated  triangular 
distribution  ranging  over  the  values  average  target  course  +  maximum  variation  in  tar¬ 
get  course  with  the  average  target  course  in  the  center  being  eight  times  as  likely  as  ei¬ 
ther  of  the  endpoints.  The  same  type  of  distribution  is  used  for  selecting  the  target's 
speed  of  advance  in  the  direction  of  the  average  target  course  on  each  leg.  These  dis¬ 
tributions  are  shown  in  Figure  2  below. 


4 


Figure  2.  Transiting  Target  Course  and  Speed  of  Advance  Distributions 


Each  leg  is  assumed  to  be  maintained  for  an  exponential  amount  of  time  as  dis¬ 
cussed  in  the  Patrolling  Target  Model  section  above.  [Ref.  3:  pp.  10-11,  20J 

flic  values  of  the  input  parameters  for  a  transiting  target  wcic  set  as  follows  to  ob 


tain  SALT  output  for  testing: 

•  Average  SOA  (knots)  12 

•  Maximum  variation  in  SOA  (knots)  4 

•  Average  target  course  (deg  N)  045 

•  Maximum  variation  in  target  course  (deg)  30 

•  Average  time  between  course  changes  (hours)  6 

C.  FLEEING  DAI  UM  TARGET  MODEL 


‘1  he  model  for  fleeing  datum  target  motion  assumes  that  the  taiget  is  clearing  datum 
from  an  area  where  it  suspects  it  has  been  detected.  To  this  end,  the  target  selects  a 
course  and  speed  from  the  specified  distributions  and  maintains  that  course  and  speed 
throughout  the  encounter.  The  target  couise  is  a  draw  from  a  Uniform  (0,  360)  distrib¬ 
ution.  The  target  speed  is  selected  from  a  truncated  triangular  distribution  with  four 
input  parameters.  The  parameters  to  be  input  are  average  speed,  minimum  speed, 
maximum  speed,  and  average  speed  likelihood  factor.  The  fust  three  parameters  are 
self-explanatory,  and  the  average  speed  likelihood  factor  determines  how  many  times 
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more  likely  the  average  speed  is  over  the  endpoint  speeds.  These  distiibutions  are  shown 
in  Figure  3  on  page  6. 


Figure  3.  Fleeing  Datum  Target  Course  anti  Speed  Distributions 

Since  the  target  maintains  course  and  speed  the  entire  time,  there  is  no  requirement  for 
a  mean  time  between  course  changes  parameter.  [Ref.  3:  pp.  10-11,  20] 

The  values  of  the  input  parameters  for  a  fleeing  datum  target  wcic  set  as  follows  to 


obtain  SALT  output  for  testing: 

•  Average  target  speed  (knots)  12 

•  Maximum  speed  (knots)  15 

•  Minimum  speed  (knots)  9 

•  Average  speed  likelihood  factor  8 


This  completes  the  discussion  of  the  three  types  of  target  motion  models  available 
for  user  selection  in  the  SALT  algotithm.  'I  he  next  chapter  will  discuss  the  development 
and  evaluation  of  sonobuoy  patterns  to  be  deployed  against  each  specific  tv  pc  of  target 
motion. 
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III.  SALT  SONOBUOY  PATTERNS 


The  sonobuoy  patterns  recommended  by  SALT  for  each  type  of  target  motion 
model  are  presented  here  in  brief  overview  form.  For  each  type  of  target  motion,  both 
the  pattern  development  and  evaluation  of  probability  of  detection  of  that  pattern  are 
discussed.  The  overall  objective  of  the  sonobuoy  patterns  developed  by  the  initial  search 
planning  phase  of  SALT  is  to  maximize  the  probability  of  detecting  the  target. 

A.  PATROLLING  TARGET  PATTERNS 

A  target  which  is  patrolling  in  an  area  as  discussed  in  the  previous  chapter  will  have 
an  average  velocity  of  zero  and  will  tend  to  remain  in  the  vicinity  of  the  initial  Search 
Probability  Area  (SPA).  The  optimal  search  density  along  each  of  the  axes  of  the 
bivariate  normal  distribution  would  be  parabolic.  The  optimization  of  search  effort 
would  require  solving  a  very  complicated  nonlinear  function  and  the  density  of  search 
effort  required  would  be  very  difficult  to  obtain  with  sonobuoys.  In  order  to  overcome 
the  difficulties  associated  with  attempting  to  optimize  the  search  effort,  the  SALT  algo¬ 
rithm  uses  an  approximation  to  the  optimal  search  effort.  The  size  of  a  rectangle  which 
most  closely  approximates  the  bivariate  normal  distribution  of  target  location  is  first 
determined.  The  available  search  effort  is  then  distributed  uniformly  over  this  rectangle. 
This  process  is  described  in  more  detail  below.  [Ref.  2:  p  36-37] 

As  the  search  for  the  target  progresses  without  any  detections,  the  ■  iiptical  area  of 
unccriaint>  of  the  initial  SPA  continues  to  expand.  To  capture  a  static  value  for  the  size 
of  the  SPA  to  use  in  detei  mining  and  evaluating  the  sonobuoy  dcplovment  pattern,  the 
initial  SPA  is  first  expanded  to  co:  respond  to  the  time  which  is  half  way  through  the 
input  search  duration.  This  expansion  gives  larger  \aiues  for  c,  and  n}  which  repicscnt 
the  average  growth  in  unccrtaintv  in  the  targets  position  over  the  duration  of  the  search. 
[Ref.  2:  p.  20-22[ 

The  bivariate  normal  distribution  dcteimincd  by  the  SPA,  expanded  as  described 
above,  can  be  expressed  as  the  product  of  two  independent  univariate  normal  distrib¬ 
utions  defined  on  the  principal  axes  of  the  SPA.  Let  a,  and  a,  be  the  standaul  deviation 
of  the  target  distiibution  along  the  first  and  seco-  1  axes.  The  lengths  of  the  sides  of  the 
rectangle  parallel  to  these  axes  arc  2 Kcr,  and  2Krr,  respectively  for  some  constant  K. 
These  values  determine  a  famih  of  rectangles  with  the  parameter  K  detcimining  the  si/c 
of  the  rectangle  as  illustrated  in  1  iguic  4  below.  [Ref.  2:  pp.  37-3S| 
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Figuie  4.  Family  of  Rectangle  Plans 


To  determine  the  probability  of  detection  for  a  rectangular  pattern  of  a  given  size, 
it  is  necessary  to  multiply  the  probability  that  the  target  is  located  in  the  lcctangle  by 
the  conditional  probability  that  the  target  is  detected  given  that  it  is  in  the  rectangle. 
'1  his  conditional  probability  is  approximated  by  an  exponential  detection  function  which 
is  based  on  the  pessimistic  assumptions  of  a  random  search  [Ref.  4:  p.  28|.  This  idea 
of  removing  the  conditioning  is  used  by  the  initial  search  planning  phase  of  SALT  to 
dctciminc  that  value  of  K  for  which  the  piobability  of  detection  calculated  by  the  fol¬ 
lowing  equation  is  maximized. 


Pd(K)  =  [<I>(K)  -  <I>( -  exp(-^)]  (3.1) 

4I\  CT|(72 


where: 

n  =  the  number  of  sonobuoys 

W  =  the  sweep  width  of  each  sonobuov  in  nautical  miles 
V  =  average  speed  of  target  as  input  by  the  user 
T  =  search  duration  in  hours 


This  maximum  value  is  output  by  SALT  as  the  probability  of  detection  for  the  recom¬ 
mended  sonobuoy  pattern.  [Ref.  2:  pp.  39-42] 

It  should  be  noted  here  that  the  term  n\VVT  in  equation  (3.1)  can  be  thought  of  as 
the  total  area  covered  by  the  pattern  during  time  period  T  due  to  target  motion  through 
the  sonobuoy  field.  We  will  later  argue  that  a  more  accurate  probability  of  detection 
can  be  calculated  by  accounting  for  the  possibility  of  detecting  the  target  at  the  begin¬ 
ning  of  the  search. 

To  develop  this  operationally  feasible  sonobuoy  pattern,  SALT  uses  the  optimally 
sized  rectangle  calculated  above  and  distributes  the  number  of  sonobuoys  available  over 
that  rectangle  to  approximate  a  uniform  distribution  of  search  effort.  To  accomplish 
this  distribution  of  effort,  the  sonobuoys  are  positioned  in  equally  spaced  rows  with  the 
spacing  between  the  buoys  in  a  row  being  equal.  If  the  number  of  sonobuoys  is  not  di¬ 
visible  b\  the  number  of  rows,  remaining  buoys  are  allotted  to  rows  beginning  with  the 
center  row  and  working  outward.  The  sonobuoy  pattern  is  output  to  the  screen  graph¬ 
ically  along  with  the  estimate  of  probability  of  detection.  [Ref.  2;  pp.  43-50) 

An  example  of  the  graphical  sonobuoy  pattern  output  for  a  scenario  involving  16 
sonobuoys  deployed  against  a  patrolling  taigct  is  presented  in  Tigurc  5  on  page  10.  The 
inner  ellipse  represents  the  initial  SPA  and  the  cuter  ellipse  represents  the  midscarch 
duration  ellipse  used  to  perfoim  the  above  calculations.  The  optimally  sized  rectangle 
is  also  displayed  in  the  pattern. 


9 


Figure  5.  Recommended  Pattern  for  16  Sonobuoys  vs.  a,  Patrolling  Target 

C.  TRANSITING  TARGET  PATI ERNS 

A  target  which  is  assumed  to  be  transiting  has  a  preferred  direction  of  motion 
which  can  be  exploited  by  the  search  planning  algorithm  by  using  barriers  which  "walk" 
against  the  target  motion.  This  fact  makes  it  possible  for  a  line  of  sonobuoys  to  sweep 
out  a  large  rectangular  area.  The  size  of  the  area  swept  depends  on  target  speed  in  the 
direction  of  the  barrier.  The  development  and  evaluation  of  sonobuoy  patterns  for  de¬ 
ployment  against  a  transiting  target  is  similar  to  that  discussed  in  the  previous  section 
for  patrolling  targets.  The  major  difference  in  the  two  developments  is  that  the 
transiting  target  search  plans  are  constructed  in  the  relative  motion  space  of  the  mean 
velocity  of  the  target  and  then  translated  to  geographic  space.  [Ref.  2:  pp.  50-59] 
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An  example  of  the  graphical  output  of  the  recommended  deployment  pattern  for  16 
sonobuoys  against  a  transiting  target  with  input  parameters  as  discussed  in  Chapter  II 
is  shown  in  Figure  6  below.  The  ellipse  displayed  in  the  figure  represents  the  initial 
SPA.  The  ellipse  which  is  updated  to  the  midsearch  duration  for  use  in  performing  cal¬ 
culations  is  not  shown. 


Figme  6.  Recommended  Pattern  for  16  Sonobuoys  vs.  a  Transiting  Target 

C.  FLEEING  DATUM  TARGET  PATTERNS 

The  development  of  a  sonobuoy  pattern  to  be  depiojed  against  a  fleeing  datum 
target  begins  by  updating  the  initial  SPA  to  the  mid-search  time.  This  update  is  ac¬ 
complished  using  a  discrete  time  and  space  Markov  chain.  As  the  search  progresses,  the 
distribution  of  target  location  spreads  radially  outward  eventually  becoming  donut 


shaped  with  the  mean,  or  center,  of  the  SPA  becoming  a  low  probability  area.  [Ref.  2: 

p.  62] 

In  order  to  accommodate  the  distribution  for  a  fleeing  datum  target,  SALT  uses  a 
class  of  search  plans  referred  to  as  "two  rectangle  plans".  These  plans  are  defined  by  an 
inner  and  an  outer  rectangle  oriented  concentrically  such  that  the  sides  of  the  rectangles 
are  parallel  to  the  major  and  minor  axes  of  the  SPA.  The  total  search  effort  is  spread 
uniformly  over  the  annular  region  between  these  two  rectangles.  The  stated  objective 
of  SALT  is  to  select  the  sizes  of  these  rectangles  so  as  to  maximize  the  probability  of 
detecting  the  target.  [Ref.  2:  p.  62] 

The  probability  of  detection  for  these  fleeing  datum  sonobuoy  patterns  is  computed 
in  a  manner  similar  to  the  calculation  discussed  above  for  a  patrolling  target.  The 
probability  that  the  target  is  contained  in  the  two  rectangle  annulus  is  multiplied  by  the 
conditional  probability  of  detecting  the  target  given  that  it  is  contained  in  this  region. 
This  conditional  probability  is  determined  by  an  exponential  detection  function  as  was 
done  in  the  patrolling  target  case.  [Ref.  2:  p.  65] 

An  example  of  this  type  of  pattern  for  16  sonobuoys  as  output  by  SALT  is  shown 
in  Figure  7  on  page  13.  The  sonobuoys  in  this  pattern  are  located  so  far  from  the  initial 
SPA  due  to  the  adjustment  of  an  input  parameter  to  the  SALT  scenario  known  as 
time-late  which  is  the  elapsed  time  between  when  the  SPA  was  obtained  and  the  time 
at  which  the  search  begins.  In  this  scenario,  time-late  has  been  set  to  four  hours. 

In  performing  the  above  calculations  of  probability  of  detection  for  the  sonobuoy 
patterns  to  be  deployed  against  any  of  the  three  types  of  assumed  target  motion,  SALT 
assumes  that  the  sonobuoys  are  all  dropped  in  place  instantly  at  the  beginning  of  the 
search.  This  assumption  does  not  adversely  aflect  the  results  of  either  the  patrolling 
target  or  'he  transiting  target.  In  the  case  of  the  patrolling  target,  the  target  tends  to 
remain  in  the  search  area  for  the  duration  of  the  search  thus  being  exposed  to  the  entire 
pattern  as  if  it  were  laid  instantaneously.  The  transiting  target  sonobuoy  patterns  are 
generally  short,  straight  line  patterns  that  do  not  take  a  long  time  to  lay.  This  fact 
combined  with  the  time  required  for  the  target  to  intercept  the  pattern  tend  to  make  this 
assumption  plausible.  1  he  sonobuoy  patterns  deployed  against  the  fleeing  datum  target, 
however,  are  usually  laid  a  large  distance  from  the  initial  target  SPA  with  large  spacing 
between  the  sonobuoys.  Due  to  the  length  of  time  required  to  lay  these  patterns  and  the 
higher  speeds  associated  with  fleeing  datum  targets,  the  assumption  of  having  the 
sonobuoys  in  place  initially  is  faulty.  During  the  data  collection  runs  of  the  fleeing  da¬ 
tum  target  scenarios,  the  target  was  sometimes  seen  to  be  beyond  the  reach  of  the  pat- 
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rigure  7.  Recommended  Pattern  for  16  Sonobuoys  vs.  a  rieeing  Datum  Target 

tern  by  the  time  the  critical  portion  of  the  pattern  was  put  in  place.  Since  the  goal  of 
this  thesis  is  only  to  test  the  validity  of  the  probability  of  detection  as  output  by  SALT, 
this  same  assumption  has  been  built  into  the  simulations  discussed  in  the  following 
chapter. 
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IV.  SIMULATION  DOCUMENTATION 


Two  examples  of  the  FORTRAN  simulation  routines  used  to  test  the  probability 
of  detection  of  the  sonobuoy  patterns  recommended  by  SALT  as  well  as  to  determine 
if  better  patterns  could  be  developed  are  presented  in  Appendix  B.  Seven  FORTRAN 
programs  are  required  to  accomplish  the  desired  testing.  For  each  target  motion  model 
there  are  two  separate  programs.  The  first  of  these  programs  is  written  to  be  used  with 
Definite  Range  Law  sonobuoys,  and  the  second  program  simulates  sonobuoys  with  de¬ 
tection  ranges  that  vary  according  to  a  Lambda-Sigma  Jump  model.  (SALT  will  handle 
either  case.)  The  remaining  program  is  a  generic  detection  subroutine  (DTCT)  which  is 
used  by  all  of  the  simulations.  Each  of  these  programs  is  discussed  below'.  All  of  the 
simulation  routines  are  developed  using  the  SALT  assumption  that  each  sonobuoy  is 
placed  in  its  exact  location  at  the  beginning  of  the  search. 

A  series  of  calibration  runs  has  been  executed  for  each  of  the  target  motion  simu¬ 
lation  programs  discussed  below  to  ensure  proper  operation.  This  testing  was  accom¬ 
plished  by  having  the  program  write  the  endpoints  of  each  track  segment  to  an  output 
file.  This  output  was  then  graphed  and  compared  to  the  respective  target  motion  model. 
Each  simulation  showed  target  motion  consistent  with  the  assumptions  for  that  type  of 
target. 

A.  PATROLLING  TARGET  SIMULATIO  S 

The  simulation  routine  for  a  patrolling  targ  with  Definite  Range  Law  sonobuoys 
is  named  "PATROL".  The  first  step  in  the  program  is  to  read  in  the  input  parameters 
from  a  data  file.  The  eight  parameters  needed  by  the  simulation  are  as  follows:  ox  and 
<jy  ,  the  standard  deviations  of  the  normal  distributions  required  to  simulate  the  initial 
Search  Probability  Area  (SPA)  and  its  containment  percent,  search  duration  in  hours, 
target's  course  change  rate  in  inverse  hours,  detection  range  of  the  sonobuoys  in  nautical 
miles,  an  integer  seed  for  use  in  the  random  number  generator,  average  target  speed  in 
knots,  and  finally,  the  Cartesian  coordinates  of  the  location  of  each  of  the  sonobuoys  in 
the  pattern. 

The  target  parameters  are  simulated  as  follows.  The  initial  target  position  is  deter¬ 
mined  by  random  draws  from  two  independent  normal  distributions  both  with  mean 
zero  and  standard  deviations  of  ax  and  or.  The  target  course,  in  radians,  is  selected  from 
a  Uniform  (0,  2n)  distribution.  The  tiiangular  distribution  required  for  the  target  speed 
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is  simulated  by  converting  a  Uniform  (0,  1)  random  number  into  a  speed  using  the  in¬ 
verse  CDF  technique.  The  speeds  generated  by  this  technique  are  displayed  as  a 
histogram  in  Figure  8  on  page  15.  It  can  be  seen  from  this  Figure  that  the  speeds  for 
patrolling  targets  in  this  simulation  are  in  fact  representative  of  a  triangular  distribution. 


'I  he  time  that  the  target  will  spend  on  the  chosen  course  at  the  chosen  speed  is  deter¬ 
mined  by  a  random  draw  from  an  exponential  distribution  with  mean 
(course  change  rate)  1  .  These  three  parameters,  course,  speed,  and  tunc,  determine  a 
line  segment  which  repiescnts  a  leg  of  the  target  track.  The  final  target  position  for  this 
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leg  of  the  track  is  calculated  using  the  initial  position,  course,  speed,  and  time  on  leg 
determined  as  discussed  above. 

The  coordinates  of  a  sonobuoy  in  the  pattern,  the  detection  range  of  that  sonobuoy, 
and  the  endpoints  of  the  current  leg  of  the  target  track  are  used  in  a  call  to  the  DTCT 
subroutine  to  determine  whether  or  not  a  detection  has  occurred.  This  process  is  re¬ 
peated  until  all  sonobuoys  in  the  pattern  under  investigation  have  been  checked  for 
possible  detections.  If  a  detection  has  occurred,  the  time  of  the  detection  is  calculated, 
and  the  number  of  detections  is  incremented  by  one.  The  detection  time  of  that 
sonobuoy  which  detects  the  target  first  in  the  time  sequence  is  written  to  an  output  file 
for  use  in  further  analysis,  and  this  run  of  the  simulation  is  terminated.  If  no  detections 
take  place  on  this  leg,  the  final  position  of  this  leg  is  reset  to  be  the  initial  position  on 
the  next  leg  and  the  time  that  the  search  has  been  underway  is  incremented.  If  this  time 
is  less  than  the  search  duration  input  by  the  user,  a  new  leg  and  final  position  are  se¬ 
lected  as  above  and  the  process  is  repeated. 

The  program  which  simulates  a  patrolling  target  in  a  sonobuoy  field  in  which  de¬ 
tection  range  is  fluctuating  according  to  a  Lambda-Sigma  Jump  model  is  entitled 
"LSPAT".  The  input  parameters  for  this  program  are  the  same  as  those  discussed  above 
with  one  exception.  Instead  of  a  constant  value  of  sonobuoy  detection  range  being  in¬ 
put,  a  Figure  of  Merit  (FOM)  is  input  by  the  user  in  dB. 

The  target  parameters,  course,  speed,  and  time  on  leg,  are  calculated  as  described 
above.  However,  a  new  environmental  parameter  is  introduced  into  this  model.  This 
parameter  is  the  acoustic  fluctuation  in  dB  which  is  represented  by  a  random  draw  from 
a  Normal  distribution  with  mean  zero  and  standard  deviation,  sigma  =  8  dB.  This 
acoustic  fluctuation  remains  constant  for  a  period  of  time  determined  from  an  expo¬ 
nential  distribution  with  mean,  lambda  =  1  hour. 

1  he  acoustic  fluctuations  simulated  as  discussed  in  the  previous  paragraph  are  con¬ 
sidered  to  be  ambient  acoustic  fluctuations  which  apply  uniformly  to  all  sonobuoys  in 
the  pattern.  The  latest  acoustic  fluctuation  is  added  to  the  input  FOM  to  determine  the 
current  FOM.  The  current  FOM  is  used  to  determine  the  detection  range  of  the 
sonobuoys  by  linear  interpolation  between  the  two  nearest  points  on  the  transmission 
loss  curve  which  is  input  as  a  separate  data  file.  The  detection  range  calculated  in  this 
way  applies  to  all  of  the  sonobuoys  in  the  pattern. 

Calls  to  the  DTCT  subroutine  for  each  sonobuoy  to  determine  detection/non- 
detcction  are  treated  as  discussed  abose.  The  only  difference  is  that  the  current  segment 
of  the  target  track  has  its  final  position  cither  at  the  end  of  the  acoustic  fluctuation  time 
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or  at  the  end  of  the  time  on  the  current  leg  whichever  occurs  next  in  the  time  sequence. 
If  the  completion  of  acoustic  fluctuation  time  occurs  first,  then  time  on  leg  is  decre¬ 
mented  by  the  acoustic  fluctuation  time  and  total  elapsed  search  time  is  incremented  by 
the  same  amount  of  time.  A  new  acoustic  fluctuation  and  time  for  the  fluctuation  are 
then  drawn  .  If  the  completion  of  the  time  on  leg  occurs  first,  then  the  acoustic  fluctu¬ 
ation  time  is  decremented  by  leg  time  and  total  time  is  incremented  by  it.  New  values 
of  course,  speed,  and  time  on  leg  are  then  computed. 

If  detection  occurs,  the  number  of  detections  is  incremented  by  one,  the  time  of  the 
detection  is  written  to  an  output  file,  and  this  run  of  the  simulation  is  finished.  If  no 
detections  have  occurred  and  the  total  elapsed  search  time  is  less  than  the  search  dura¬ 
tion  input  by  the  user,  the  process  continues. 

B.  TRANSITING  TARGET  SIMULATION 

The  program  which  simulates  a  transiting  target  versus  a  field  of  Definite  Range 
Law  sonobuoys  functions  in  a  manner  similar  to  that  described  in  the  above  section  on 
Patrolling  Target  Simulations.  The  input  parameters  for  this  program  are  the  same  as 
those  discussed  above  except  that  the  minimum  and  maximum  course  are  input  in 
„  adians  and  the  minimum  and  maximum  speed  are  input  in  knots. 

The  transiting  target  parameters  are  simulated  as  follows.  The  initial  position  is 
determined  the  same  way  as  described  under  the  Patrolling  Target  Simulations  section 
above.  Both  target  course  and  speed  of  advance  (SOA)  are  simulated  from  triangular 
distributions  using  the  inverse  CDF  technique.  These  values  of  course  and  SOA  are  then 
converted  to  target  speed  by  dividing  the  SOA  by  the  cosine  of  the  angle  between  aver¬ 
age  target  course  and  the  target  course  determined  for  this  track  segment.  These  trian¬ 
gular  distributions  have  been  verified  by  histogram  as  in  the  patrolling  target  case,  but 
the  figure  is  not  reproduced  here.  The  remainder  of  the  simulation  proceeds  exactly  as 
described  above  for  patrolling  targets. 

The  program  which  models  the  Lambda-Sigma  Jump  sonobuoys  generalizes  from 
the  transiting  target  Definite  Range  Law  sonobuoy  program  in  the  same  manner  as  for 
patrolling  targets. 

C.  FLEEING  DATUM  TARGET  SIMULATION 

The  next  program  to  be  discussed  is  the  one  used  to  simulate  fleeing  datum  targets 
and  Definite  Range  Law  sonobuoys.  The  input  parameters  are  the  same  as  those  dis¬ 
cussed  in  the  previous  sections.  The  target  course  is  drawn  from  a  Uniform  (0,  2n)  dis¬ 
tribution  and  target  speed  is  drawn  from  a  triangular  distribution  determined  by  the 
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input  parameters  of  minimum  and  maximum  target  speed.  The  course  and  speed  need 
only  be  drawn  one  time  for  each  simulation  run  since  the  target  is  assumed  to  maintain 
that  course  and  speed  throughout  the  engagement.  The  initial  position  is  aged  by  the 
time-late  value  of  four  hours  using  the  selected  course  and  speed.  The  detection/non- 
detection  event  is  then  determined  as  described  in  the  Patrolling  Target  Simulations 
section.  This  simulation  routine  generalizes  to  the  Lambda-Sigma  Jump  model  for 
sonobuoy  detection  range  in  the  same  manner  as  previously  discussed. 

D.  DETECTION  SUBROUTINE 

The  detection  subroutine  entitled  "DTCT"  is  used  by  all  six  of  the  main  programs 
to  determine  whether  or  not  a  target  detection  occurs.  If  a  detection  does  take  place, 
the  subroutine  returns  the  distance  from  the  beginning  of  the  current  segment  of  the 
target  track  to  the  point  where  the  earliest  detection  occurs. 

The  inputs  to  the  DTCT  subroutine  aie  the  detection  range  of  the  sonobuoys.  the 
cartesian  coordinates  of  the  endpoints  of  the  current  track  segment,  and  the  location  of 
the  current  sonobuoy  being  tested  for  detection  criterion.  These  seven  parameters  are 
passed  into  the  subroutine  at  each  call  to  it. 

The  first  step  of  this  subroutine  is  to  calculate  the  perpendicular  distance  between 
the  line  determined  by  the  target's  track  segment  and  the  location  of  the  sonobuoy.  If 
this  value  is  greater  than  the  detection  range  of  the  sonobuoy,  no  detection  takes  place. 
However,  if  this  perpendicular  distance  is  less  than  the  detection  range,  there  arc  three 
cases  to  consider  to  determine  wnethcr  ^r  not  detection  occurs.  If  the  track  segment 
contains  the  point  on  the  line  that  is  closest  to  the  buoy  center,  then  a  detection  occurs. 
If  the  track  segment  does  not  contain  that  point,  but  either  endpoint  of  the  track  seg¬ 
ment  falls  within  the  range  of  detection  of  the  buoy  a  detection  will  also  take  place.  If 
neither  of  these  conditions  holds,  then  there  is  no  detection  for  this  track  segment- 
sonobuoy  combination.  Each  of  these  four  cases  is  illustrated  in  Tigure  9  on  page  19. 

In  order  to  ensure  the  proper  operation  of  this  subroutine,  examples  of  the  four 
cases  described  above  were  developed  using  simple  graphical  techniques.  These  exam¬ 
ples  were  then  input  into  the  subroutine  for  testing.  The  outputs  from  the  subroutine 
confirm  that  it  performs  its  function  properly  as  described  here. 
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Figure  9. 


Sonobuoy  Detection/Non-detection  Possibilities 
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V.  SIMULATION  RESULTS  OF  SALT  SONOBUOY  PATTERNS 


The  sonobuoy  patterns  generated  by  the  SALT  algorithm  for  each  specific  type  of 
target  motion  were  tested  using  the  respective  simulation  routine  discussed  in  the  previ¬ 
ous  chapter  to  determine  if  the  probability  of  detection  for  that  pattern  as  output  by 
SALT  is  accurate.  The  simulation  results  are  output  as  a  number  of  detections  which 
can  be  easily  converted  to  a  probability  of  detection  if  the  number  of  detections  is  di¬ 
vided  by  the  number  of  runs  of  the  simulation  which  is  10,000  runs  in  each  case.  Each 
proportion,  p,  which  represents  the  probability  of  detection  can  be  expanded  to  a  two- 
sided  (1  —  a)%  Confidence  Interval  by  use  of  the  following  formula: 


where  Z,_i  represents  the  proper  percentage  point  of  a  Standard  Normal  cumulative 
distribution. 

The  SALT  algorithm  was  run  three  times  for  each  type  of  target  motion  using 
sonobuoy  s  that  have  a  Definite  Range  Law  detection  function.  The  parameters  within 
the  SALT  program  which  determine  the  acoustic  fluctuations  were  set  to  their  minimum 
values  to  accomplish  this  (see  Appendix  A).  For  each  run  of  a  specific  target  motion 
type,  all  inputs  remain  the  same  with  the  exception  of  the  number  of  sonobuoys.  The 
number  of  sonobuoys  used  for  each  run  are  8,  16,  and  24.  These  numbers  of  sonobuoys 
allowed  for  testing  over  a  range  of  values  of  probability  of  detection.  The  same  nine 
runs  of  the  SALT  algorithm  were  repeated  to  obtain  data  for  sonobuoys  which  behave 
according  to  a  Lambda-Sigma  Jump  model.  The  parameters  for  this  phase  of  testing 
are  an  ambient  acoustic  fluctuation  rate  of  one  per  hour  and  standard  deviation  of  am¬ 
bient  acoustic  fluctuations  of  eight  dB.  The  probability  of  detection  as  output  by  SALT 
will  now  be  compared  to  the  simulation  results  for  each  type  of  target  motion. 

A.  PATROLLING  TARGET 

The  probability  of  detection  for  the  Definite  Range  Law  sonobuoy  patterns  recom¬ 
mended  by  SALT  for  a  patrolling  target  are  displayed  in  Table  1  on  page  21.  Both  the 
SALT  calculated  probability  of  detection  and  the  simulation  routine  output  probability 
of  detection  are  listed  for  the  number  of  sonobuoys  used  in  each  pattern.  A  95%  Con- 


fidence  Interval,  calculated  as  described  above,  is  also  displayed  for  each  simulation 
probability  of  detection. 


Table  1.  PATROLLING  TARGET  COMPARISONS  (DEFINITE  RANGE) 


NUMBER 

SALT 

SIMULATION 

95%  CONFIDENCE 

OF  BUOYS 

Pd 

Pd 

INTERVAL 

8 

0.49 

0.8848 

(0.8785.  0.S911) 

16 

0.67 

0.9402 

(0.9378,  0.9426) 

24 

0.76 

0.9773 

(0.9758,  0.9788) 

It  can  be  seen  from  this  table  that  the  SALT  algorithm  tends  to  understate  the 
probability  of  detection  for  a  target  in  this  case.  The  low  probabintics  of  detection 
output  by  SALT  have  the  potential  to  be  very  misleading.  For  example,  the  SALT  re¬ 
commended  pattern  for  eight  sonobuoys  lists  as  its  probability  of  detection  a  value  of 
0.49.  This  low  probability  of  detection  may  encourage  a  P-3C  crew  to  opt  instead  for  a 
sonobuoy  pattern  with  more  buoys  which  would  be  a  waste  of  resources  considering  the 
actual  probability  of  detection  for  the  target  using  an  eight  buoy  pattern  is  0.88. 

The  Technical  Documentation  of  SALT  claims  that  the  exponential  detection  func¬ 
tion  used  in  equation  (3.1)  on  page  9  gives  a  lower  bound  on  the  probability  of  detection 
that  can  be  achieved  by  a  systematic  uniform  search  of  a  region  [Ref.  2:  p.  26].  Due  to 
this  fact,  one  would  expect  to  get  resulting  probabilities  of  detection  that  are  a  little 
higher  than  the  SALT  output  values.  However,  the  disparity  in  the  values  given  here  is 
too  high  to  be  so  easily  explained  away.  In  order  to  gain  insight  into  why  the  simulation 
probability  of  detection  is  so  different  from  the  SALT  probability  of  detection,  the  de¬ 
tection  times  as  output  by  the  simulation  program  are  graphed  into  a  Cumulative  Dis¬ 
tribution  Function  in  Figure  10  on  page  22. 
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Figure  10.  Cumulative  Distribution  of  Detection  limes 

It  can  be  seen  from  this  CDF  that  a  large  percentage  of  the  targets  which  get  de¬ 
tected  arc  initial  detections.  This  is  not  too  surprising  considering  the  tight  packing  of 
the  sonobuoys  into  the  86%  containment  ellipse.  The  sonobuoy  detection  circles  are 
shown  plotted  over  this  86%  containment  ellipse  in  Figure  1 1  on  page  23,  The  ellipse 
shown  in  this  figure  has  a  semi-major  axis  of  60  nautical  miles  and  a  semi-minor  axis  of 
30  nautical  miles.  The  detection  range  circles  are  18.2  nautical  miles  in  radius  and  are 
placed  in  the  locations  specified  by  SALT. 
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Figure  11.  Detection  Ranges  and  86%  Containment  Ellipse  for  an  8  Sonobuoy 
Pattern 

To  account  for  this  initial  probability  of  detection  in  the  patrolling  target  sonobuoy 
patterns,  equation  (3.1)  on  page  9  can  be  modified  slightly  to  read  as  follows: 

Pd(K)  =  [<D(K)  -  cb(  -K)]2[l  -  exp(  -  .nWVl+-P..gRl )]  (5.2) 

4k  o,o2 

where: 

R  =  the  detection  range  of  the  sonobuoys  in  nautical  miles 

and  all  other  variables  are  defined  as  they  are  in  equation  (3.1). 

The  above  modification  to  the  calculation  of  the  probability  of  detection  takes  into 
account  the  detection  area  of  each  sonobuoy  at  the  initiation  of  search.  The  value  of 
the  product  nWVT  is  clearly  time  dependent.  This  means  that  in  equation  (3.1)  there 
is  zero  probability  of  an  initial  detection.  The  detection  times  shown  in  Figure  10  on 
page  22  clearly  show  the  need  to  account  for  this  initial  pi  obability  of  detection. 


Using  equation  (5.2)  to  calculate  the  probability  of  detection  for  the  three  Definite 
Range  Law  sonobuoy  patterns  specified  in  Table  1  on  page  21  results  in  probabilities 
of  detection  of  0.74,  0.89,  0.93  for  8,  16,  and  24  sonobuoy  patterns  respectively.  While 
these  values  do  not  exactly  reflect  the  results  of  the  simulation  runs,  they  do  provide  the 
decision  maker  with  a  much  more  realistic  assessment  of  the  probability  of  detection  for 
the  recommended  pattern. 

The  Lambda-Sigma  Jump  model  comparisons  for  the  similar  patterns  are  given  in 
Table  2.  The  results  for  the  8  and  16  sonobuoy  patterns  show  the  same  phenomenon 
discussed  above  for  Definite  Range  Law  sonobuoys.  In  these  cases,  equation  (5.2)  again 
provides  a  more  realistic  value  of  probability  of  detection. 


Table  2.  PATROLLING  TARGET  COMPARISONS  (LAMBDA-SIGMA) 


NUMBER 

OF  BUOYS 

SALT 

Pd 

SIMULATION 

PD 

95%  CONFIDENCE 
INTERVAL 

8 

0.59 

0.8664 

(0.8630.  0.8698) 

16 

0.76 

0.9028 

(0.8994.  0.9062) 

24 

0.84 

0.8120 

(0.8081.  0.8159) 

The  SALT  output  for  24  sonobuoys  in  this  scenario  was  a  pattern  which  was  ro- 
ticeably  different  from  the  five  other  sonobuoy  patterns  recommended  for  deployment 
against  a  patrolling  target.  The  SALT  generated  pattern  is  a  BRUSHTAC  pattern  where 
all  of  the  buoys  are  placed  in  a  narrowly  spaced  line  along  the  major  axis  of  the  initial 
SPA.  In  a  telephone  conversation  with  Lawrence  I).  Stone,  the  president  of  Metron, 
Inc.,  it  has  been  determined  that  an  error  in  the  SALT  code  resulted  in  an  optimal  target 
depth  being  selected  which  would  give  the  sonobuoys  a  convergence  zone  detection. 
This  selection  was  made  by  the  algorithm  in  spite  of  the  fact  that  the  probability  of  the 
target  being  at  that  depth  was  input  as  zero.  Due  to  this  unfortunate  fact,  the  results 
displayed  in  Table  2  for  a  24  sonobuoy  pattern  have  no  meaning. 

B.  TRANSITING  TARGET 

'Ihc  comparison  of  the  SALT  output  probability  of  detection  for  transiting  target 
sonobuoy  patterns  to  the  probability  of  detection  for  the  same  patterns  as  output  by  the 
simulation  routine  is  presented  in  Table  3  on  page  25. 
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Table  3.  TRANSITING  TARGET  COMPARISONS  (DEFINITE  RANGE) 


NUMBER 

OF  BUOYS 

SALT 

Pd 

SIMULATION 

Pd 

95%  CONFIDENCE 
INTERVAL 

8 

0.72 

0.8857 

(0.8825,  0.8889) 

16 

0.81 

0.9085 

(0.9056,  0.9114) 

24 

0.82 

0.9156 

(0.9128,0.9184) 

The  probabilities  of  detection  displayed  above  do  not  exhibit  the  same  magnitude 
of  difference  that  those  of  the  patrolling  target  exhibit  in  the  previous  section  as  seen  in 
Table  1  on  page  21.  However,  there  is  still  a  slight  difference  in  these  probabilities  of 
detection  that  requires  consideration.  The  CDF  plot  for  transiting  target  detection  times 
reveals  that  only  about  35%  of  the  detections  obtained  are  due  to  initial  detection.  The 
reason  for  this  smaller  number  of  initial  detections  is  that  the  transiting  target  search 
plans  are  constructed  in  the  relative  motion  space  of  the  mean  velocity  of  the  target  as 
discussed  in  the  Transiting  Target  Patterns  section  of  Chapter  III.  The  length  of  the  side 
of  the  rectangular  area  swept  out  which  is  parallel  to  the  mean  velocity  vector  is  deter¬ 
mined  by  the  following  equation:  [Ref.  2:  p.  52] 

Length  =  VT  +  2R  (5.3) 


where: 

V  -  average  SOA  of  the  target  in  knots 

T  =  search  duration  in  hours 

R  =  the  detection  range  of  the  sonobuoys  in  nautical  miles 

When  this  relative  motion  space  is  translated  into  geographic  space,  only  the  2R 
portion  of  the  length  as  determined  by  equation  (5.3)  has  the  possibility  for  initial  de¬ 
tections.  In  order  to  account  for  these  initial  detections  in  the  calculation  of  probability 
of  detection,  this  fraction  of  the  total  length  should  be  used. 

Due  to  the  same  depth  selection  error  discussed  at  the  end  of  the  previous  section, 
the  SALT  output  for  transiting  targets  using  the  Lambda-Sigma  Jump  model  is  only 
valid  for  the  eight  sonobuoy  case.  The  SALT  output  probability  of  detection  for  this 
scenario  is  0.S1  while  the  simulation  output  for  the  identical  pattern  is  0.7996.  The  95"  o 


Confidence  Interval  about  this  value  is  (0.7956,  0.8036).  In  this  more  realistic  case,  the 
initial  probability  of  detection  does  not  seem  to  be  an  important  factor. 

C.  FLEEING  DATUM  TARGET 

The  results  of  the  simulation  analysis  of  Definite  Range  Law  sonobuoy  patterns 
deployed  against  a  fleeing  datum  target  are  shown  in  Table  4. 


Table  4.  FLEEING  DATUM  TARGET  COMPARISONS  (DEFINITE  RANGE) 


NUMBER 

OF  BUOYS 

SALT 

Pd 

SIMULATION 

Pd 

95%  CONFIDENCE 
INTERVAL 

8 

0.47 

0.4829 

(0.4779,  0.4879) 

16 

0.69 

0.9579 

(0.9559,  0.9599) 

24 

0.81 

0.9985 

(0.9981,0.9989) 

The  probability  of  detection  as  output  by  both  SALT  and  the  simulation  of  this 
scenario  appear  to  be  in  close  agreement  for  the  eight  sonobuoy  case  while  the  other  two 
patterns  show  a  highly  significant  difference.  This  difference  cannot  be  explained  by 
initial  detections  as  is  done  in  the  cases  of  patrolling  and  transiting  targets  above.  Less 
than  10%  of  the  detections  were  achieved  initially  when  these  scenarios  were  run  using 
the  simulation  routines.  An  alternate  reason  for  this  phenomenon  can  be  found  by  de¬ 
termining  the  fraction  of  the  perimeter  of  these  patterns  that  is  co\ered  by  sonobuoy 
detection  radii.  In  the  eight  sonobuoy  pattern,  48%  of  this  perimeter  is  covered  by  such 
detection  ranges.  Similar  calculations  for  the  16  sonobuoy  pattern  reveal  that  94"  o  of 
the  perimeter  is  so  covered,  and  the  percentage  of  coverage  for  the  24  sonobuoy  case  is 
100%.  This  fact  accurately  reflects  that  when  the  sonobuoys  are  assumed  to  be  placed 
in  the  water  immediately  upon  commencement  of  search,  it  is  highly  unlikely  that  a 
target  will  avoid  having  to  pass  through  the  annular  region  where  the  sonobuoys  are 
placed. 

As  discussed  in  the  Fleeing  Datum  Target  Patterns  section  of  Chapter  III,  the 
probability  of  detection  in  these  scenarios  is  calculated  within  the  SALT  algorithm  by 
muitiplving  the  probability  of  detecting  the  target  given  that  it  is  contained  in  the  two 
rectangle  annular  region  by  the  probability  that  the  target  is  contained  in  this  region. 
By  performing  this  simplified  calculation,  SALT  fails  to  take  into  account  the  fact  that 
a  target  which  is  inside  of  the  inner  rectangle  at  the  commencement  of  search  must  pass 
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through  this  annular  region  to  evade  detection.  The  SALT  algorithm  must  be  modified 
in  some  way  to  account  for  this  possibility  in  order  to  more  realistically  evaluate  the 
probability  of  detection  for  fleeing  datum  targets. 

The  results  for  the  same  target  scenarios  versus  Lambda-Sigma  Jump  model 
sonobuoys  are  presented  in  Table  5.  These  results  show  a  disparity  similar  to  that  dis¬ 
cussed  above  for  Definite  Range  Law  sonobuoys.  Thus  the  corrective  action  specified 
above  is  valid  for  both  types  of  sonobuoy  detection  models. 


Table  5.  FLEEING  DATUM  TARGET  COMPARISONS  (LAMBDA-SIGMA) 


NUMBER 

OF  BUOYS 

SALT 

Pd 

SIMULATION 

Pd 

95%  CONFIDENCE 
INTERVAL 

8 

0.59 

0.8623 

(0.8589.  0.8657) 

16 

0.80 

0.9138 

(0.9110,  0.9166) 

24 

0.88 

0.9413 

(0.9389,  0.9437) 

It  can  be  seen  from  these  analyses  that  in  all  cases  of  target  motion  models  the 
SALT  code  requires  modification  to  be  able  to  more  accurately  reflect  the  probability 
of  detection  so  that  sound  decisions  can  be  made  by  the  user  of  this  tactical  decision  aid. 
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VI.  ALTERNATIVE  SONOBUOY  PATTERN  DEVELOPMENT 

The  second  phase  of  this  research  effort  involves  attempting  to  determine  whether 
or  not  the  sonobuoy  patterns  output  by  the  SALT  algorithm  can  be  improved.  The 
measure  of  effectiveness  of  a  sonobuoy  pattern  will  be  the  probability  of  detection  as 
output  by  the  same  simulation  routines  used  to  verify  the  probability  of  detection.  The 
improvements  considered  and  the  results  are  presented  for  each  target  motion  model  in 
the  following  sections. 

A.  PATROLLING  TARGET  MODEL 

As  was  shown  previously  in  Figure  1 1  on  page  23,  the  pattern  for  eight  r.onobuoys 
deployed  against  a  patrolling  target  consists  of  excessive  overlap  of  the  buoy  detection 
regions.  The  sonobuoy  spacing  in  this  pattern  is  22.3  nautical  miles  between  the  two 
rows  and  21.2  nautical  miles  between  the  buoys  within  a  row.  The  probability  of  de¬ 
tection  of  this  pattern  is  0.88  as  determined  from  the  simulation  program. 

The  SALT  assumption  for  deploying  a  sonobuoy  pattern  against  a  target  of  this 
type  is  that  the  search  effort  should  be  distributed  uniformly  over  the  region  in  which  the 
target  is  most  likely  to  be  located.  Under  this  assumption  it  can  plainly  be  seen  that  the 
overlap  of  the  buoy  detection  ranges  is  excessive.  By  increasing  the  spacing  between  the 
buoys  in  this  pattern,  some  gains  in  probability  of  detection  can  be  realized. 

The  best  spacing  to  use  between  the  sonobuovs  in  the  pattern  with  eight  Definite 
Range  Law  sonobuoys  has  been  determined  by  trail  and  error.  After  many  simulation 
runs  on  patterns  with  various  spacing,  it  has  been  determined  that  the  best  effort  can 
be  achieved  by  making  the  spacing  between  rows  31.2  nautical  miles,  and  setting  the 
within  row  spacing  to  36.0  nautical  miles.  The  probability  of  detection  for  this  pattern 
was  output  by  the  simulation  as  0.96.  Further  spreading  of  the  sonobuoys  results  in  too 
much  effort  being  allocated  to  outlying  areas.  The  probability  that  the  target  is  located 
in  these  areas  is  very  low. 

A  graphical  representation  of  the  detection  ranges  of  the  sonobuoys  superimposed 
on  an  86%  probability  of  containment  ellipse  is  shown  in  Figure  12  on  page  29.  The 
ellipse  shown  in  the  figure  has  a  semi-major  axis  of  60  nautical  miles  and  a  semi-minor 
axis  of  30  nautical  miles.  The  defection  range  circles  arc  IS. 2  nautical  miles  in  radius, 
and  overlap  only  a  slight  amount. 
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Figure  12.  Detection  Ranges  and  86%  Containment  Ellipse  for  an  8  Sonobuoy 
Pattern  with  36  Nautical  Mile  Spacing 


Similar  tests  were  run  using  the  Lambda-Sigma  Jump  model  sonobuoys,  but  altering 
the  spacing  docs  not  affect  the  results  as  much  since  them  arc  changes  in  sonobuoy  de¬ 
tection  range  caused  by  acoustic  fluctuations.  The  SALT  output  spacing  for  the  eight 
sonobuoy  Lambda-Sigma  Jump  pattern  is  24.7  nautical  miles  between  rows  and  23.5 
nautical  miles  between  buoys  within  each  row.  Altering  the  spacing  between  the 
sonobuoys  in  this  scenario  provides  no  significant  increases  in  the  probability  of  de¬ 
tection. 

Another  type  of  pattern  that  was  tested  using  Lambda-Sigma  Jump  model 
sonobuoys  against  a  patrolling  target  is  a  2  x  4  x  2  pattern.  This  pattern  consists  of  4 
sonobuoys  placed  along  the  major  axis  of  the  86%  containment  ellipse  with  2  sonobuoys 
on  each  side  of  this  line  located  near  the  minor  axis.  No  notable  increase  in  probability 
of  detection  is  realized  using  this  pattern. 
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Based  on  the  results  of  the  testing  described  above,  it  seems  that  the  sonobuoy  pat¬ 
terns  recommended  by  SALT  for  use  in  the  more  realistic  environmental  setting  modeled 
by  the  Lambda-Sigma  Jump  process  are  consistent  with  the  stated  objectives  of  the 
SALT  algorithm. 

B.  TRANSITING  TARGET  MODEL 

In  attempting  to  improve  the  sonobuoy  patterns  recommended  by  SALT  for  use 
against  a  transiting  target,  only  the  more  realistic  case  of  the  Lambda-Sigma  Jump 
model  sonobuoys  is  considered.  The  SALT  pattern  for  this  type  of  target  motion  model 
is  a  straight  line  of  buoys  placed  perpendicular  to  the  mean  target  course.  In  looking 
at  this  pattern  of  eight  sonobuoys  graphed  along  with  the  86%  containment  ellipse,  it 
appears  that  the  line  of  sonobuoys  is  located  too  close  to  the  center  of  the  SPA.  Moving 
the  line  of  sonobuo)  s  away  from  the  center  in  the  direction  of  the  mean  target  course 
results  in  the  probability  of  detection  going  down.  This  results  from  the  fact  that  al¬ 
though  fewer  targets  will  be  located  on  the  wrong  side  of  this  line,  more  targets  will  be 
located  further  from  the  pattern  and  will  not  reach  it  prior  to  termination  of  the  search. 
The  next  aspect  of  the  pattern  to  be  considered  is  the  spacing  between  the  sonobuoys 
in  the  line.  The  spacing  between  the  sonobuoys  in  the  SALT  generated  pattern  is  only 
12.8  nautical  miles,  even  though  the  average  detection  range  is  26.3  nautical  miles.  The 
sonobuoy  spacing  was  increased  to  20  nautical  miles.  The  simulation  run  for  this  new 
pattern  showed  no  notable  increase  in  probability  of  detection.  This  result  is  not  diffi¬ 
cult  to  understand  in  light  of  the  fact  that  as  the  spacing  between  the  sonobuo)  s  is  in¬ 
creased,  the  buoxs  near  the  ends  of  the  line  are  being  moved  to  areas  with  extremely  low 
probability  of  achieving  a  detection. 

The  next  logical  step  in  the  development  of  alternative  patterns  is  to  depart  from  the 
requirement  to  lay  the  sonobuoys  in  a  straight  line.  By  increasing  the  spacing  between 
the  sonobuo) s  to  overcome  the  excessive  overlap  and  then  bringing  the  endpoints  of  the 
line  closer  to  the  SPA,  it  seems  that  the  probability  of  detection  should  increase.  Two 
patterns  were  tested  using  this  theory. 

The  first  pattern  tested  in  which  the  sonobuoys  are  not  placed  in  a  straight  line  is  a 
semi-circular  pattern  with  spacing  between  the  sonobuoys  of  20  nautical  miles.  This 
pattern  covers  90  degrees  of  arc  on  each  side  of  the  mean  target  course  line.  The  pattern 
is  positioned  so  that  all  sonobuoys  are  located  at  a  distance  midway  between  the  lengths 
of  the  semi-major  and  semi-minor  axes  from  the  center  of  the  SPA. 
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The  second  pattern  tested  is  a  chevron  pattern.  The  sonobuoys  in  this  pattern  are 
divided  equaily  between  the  two  sides  of  the  chevron,  and  the  wedge  angle  of  the  chevron 
is  set  at  90  degrees.  The  point  of  the  chevron  is  located  along  the  mean  target  course 
line  at  a  distance  half  way  between  the  lengths  of  the  semi-major  and  semi-minor  axes. 
The  spacing  between  the  sonobuoys  was  again  set  to  20  nautical  miles. 

Neither  of  these  patterns  results  in  a  significant  increase  in  the  probability  of  de¬ 
tection  for  the  transiting  target  modeled  here.  Based  on  the  test  results  described  here, 
the  sonobuoy  patterns  generated  by  SALT  for  deployment  against  transiting  targets 
appear  to  adequately  meet  the  objectives  of  SALT. 

C.  FLEEING  DATUM  TARGET  MODEL 

The  first  sonobuoy  pattern  output  by  SALT  that  is  considered  for  improvement  is 
the  pattern  with  eight  Definite  Range  Law  sonobuoys.  The  analysis  of  the  detection 
times  output  by  the  simulation  routine  for  this  pattern  shows  only  42  initial  detections 
out  of  the  4829  detections  obtained.  This  result  leads  to  the  suspicion  that  this  pattern 
is  located  too  far  from  the  initial  SPA.  Each  side  of  the  square  pattern  generated  by 
SALT  is  located  88  nautical  miles  from  the  center  of  the  SPA.  Several  runs  of  the  sim¬ 
ulation  have  been  conducted  moving  the  sonobuoy  pattern  nearer  the  SPA  each  time 
and  scaling  it  accordingly.  These  runs  show  that  a  pattern  in  which  the  sides  are  located 
50  nautical  miles  from  the  center  increases  the  probability  of  detection  to  0.78. 

The  square  shape  of  the  sonobuoy  patterns  output  by  SALT  for  attempting  to  de¬ 
tect  a  fleeing  datum  target  naturally  suggests  that  a  better  use  of  available  search  effort 
would  be  to  arrange  the  sonobuoys  in  a  circle  so  that  the  buoys  would  be  equidistant 
from  the  SPA  center.  Such  a  pattern  was  tested  using  the  eight  sonobuoy  Lambda- 
Sigma  Jump  model  scenario.  The  sonobuoys  in  this  pattern  were  placed  82.5  nautical 
miles  from  the  center  of  the  SPA  and  arranged  so  as  to  uniformly  cover  the  perimeter 
of  the  circle.  This  distance  from  the  center  was  arrived  at  by  averaging  the  semi-major 
and  semi-minor  axes  of  the  initial  SPA  and  then  adding  a  distance  to  account  for  the 
time-late  of  the  beginning  of  the  search  which  is  four  hours.  In  this  way,  all  sonobuoys 
in  the  pattern  are  located  closer  to  the  center  than  in  the  SALT  pattern.  Each  edge  of 
the  square  SALT  pattern  is  located  at  84.6  nautical  miles  from  the  center. 

A  similar  circular  pattern  has  also  been  tested  in  which  the  sonobuoys  are  placed 
60  nautical  miles  from  the  center.  This  distance  is  based  on  the  improved  pattern  for 
Definite  Range  Law  sonobuoys  discussed  above.  This  pattern  only  increases  the  prob¬ 
ability  of  detection  by  0.03  which  is  not  significant. 
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The  above  results  of  attempting  to  improve  the  sonobuoy  patterns  recommended 
by  SALT  again  show  that  the  Definite  Range  Law  scenario  patterns  can  be  improved 
greatly  while  the  Lambda-Sigma  Jump  scenario  patterns  do  appear  to  be  consistent  with 
the  SALT  objectives.  These  results  should  not  be  weighted  too  heavily  as  they  are  ob¬ 
tained  by  remaining  within  the  constraints  of  the  SALT  assumptions,  most  notably  that 
all  of  the  sonobuoys  in  a  pattern  deployed  against  a  fleeing  datum  target  are  put  in  place 
at  the  commencement  of  the  search.  This  assumption  has  been  shown  to  be  faulty  at 
the  end  of  Chapter  III. 
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VII.  CONCLUSIONS 


It  has  been  established  at  the  beginning  of  this  thesis  that  the  goal  of  the  research 
effort  is  to  investigate  the  initial  search  planning  phase  of  the  Search  and  Localization 
Tactical  Decision  Aid.  This  investigation  consists  of  first  attempting  to  ensure  that  the 
probability  of  detection  of  various  sonobuoy  patterns  recommended  by  SALT  is  correct 
under  the  stated  assumptions  and  subsequently  of  testing  alternative  sonobuoy  patterns 
which  could  lead  to  a  better  patterns  than  those  recommended  by  SALT. 

The  probability  of  detection  output  by  SALT  for  its  own  patterns  is  significantly 
understated  for  both  the  patrolling  target  model  and  the  fleeing  datum  target  model. 
Possible  corrections  to  these  errors  have  been  presented  here.  The  probability  of  de¬ 
tection  for  the  transiting  target  model  appears  to  be  correct  especially  in  the  more  real¬ 
istic  case  of  the  Lambda-Sigma  Jump  model  for  sonobuoy  detection  ranges. 

The  sonobuoy  patterns  generated  by  SALT  for  Definite  Range  Law  sonobuoys  de¬ 
ployed  against  patrolling  and  fleeing  datum  targets  have  been  shown  to  be  inadequate 
for  accomplishing  the  goal  of  maximizing  probability  of  detection.  The  recommended 
patterns  for  all  three  motion  models  when  the  acoustic  fluctuations  which  influence 
sonobuoy  detection  range  are  modeled  by  a  Lambda-Sigma  Jump  process  are  consistent 
with  the  stated  objectives  of  SALT.  This  fact  is  especially  encouraging  since  the  patterns 
output  by  the  SALT  algorithm  are  operationally  easy  to  deploy  from  a  P-3C  aircraft. 
One  must  be  reminded  here  that  the  testing  performed  on  the  sonobuoy  patterns  output 
by  SALT  has  been  conducted  within  the  constraints  of  the  assumptions  made  in  devel¬ 
oping  the  SALT  algorithm. 

One  aspect  of  the  fleeing  datum  target  scenario  which  warrants  further  investigation 
is  the  assumption  that  all  sonobuoys  are  put  in  place  simultaneously  at  the  commence¬ 
ment  of  the  search.  Some  method  of  calculating  the  probability  of  detection  for  this  type 
of  target  needs  to  be  developed  which  allows  for  the  fact  that  a  target  may  be  outside 
of  the  sonobuoy  pattern  by  the  time  that  the  critical  portion  of  these  large  patterns  is 
put  in  place.  A  more  useful  alternative  to  altering  the  calculation  of  the  probability  of 
detection  would  be  to  develop  a  means  of  accounting  for  the  required  flight  time  in  re¬ 
commending  the  sonobuoy  pattern.  This  change  in  the  algorithm  would  open  up  the 
possibility  of  deploying  patterns  which  spiral  out  from  the  center  as  time  elapses  and  the 
target  continues  its  assumed  motion  away  from  the  SPA  center.  [Ref.  4:  pp.  115-118] 
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The  SALT  algorithm  is  rated  as  effective  in  its  initial  search  planning  phase  for  pa- 
trol'ing  targets  with  only  minor  modifications  necessary  to  improve  the  calculations  of 
probability  of  detection.  The  SALT  recommendations  from  the  initial  search  planning 
phase  for  transiting  targets  also  appear  to  be  realistic.  Further  development  is  required 
in  the  case  of  fleeing  datum  targets  as  discussed  above  in  order  to  make  SALT  ready  for 
operational  testing. 


APPENDIX  A.  SALT  INPUTS  FOR  TEST  CASES 


The  purpose  of  this  appendix  is  to  provide  a  list  of  the  inputs  to  the  SALT  algorithm 
which  were  used  to  generate  the  sonobuoy  patterns  required  for  the  testing  performed 
here.  Only  those  parameters  which  are  relavent  to  the  scenarios  are  presented.  The  re¬ 
maining  inputs  are  left  at  their  default  values.  The  target  motion  model  parameters  are 
not  presented  here  since  they  are  given  in  Chapter  II. 

The  following  list  shows  input  values  used  in  the  Definite  Range  Law  sonobuoy 
cases  for  patrolling  and  transiting  targets. 


Setup 

Buoy  acoustic  fluctuation  rate  (per  hr)  0.01 

Standard  deviation  of  buoy  acoustic  fluctuations  (dB)  0.01 

Ambient  acoustic  fluctuation  rate  (per  hour)  0.01 

Standard  deviation  of  ambient  acoustic  fluctuations  (dB)  0.01 

Search  duration  (hrs)  4.0 

Environmental  Data 

Ocean  Pacific 

Ambient  noise  level  (dB)  75 

Target  Parameters 

Depth  1  (ft)  90 

Probability  of  depth  1  1.0 

Depth  2  (ft)  300 

Probability  of  depth  2  0.0 

Frequency  1  (hz)  50 

Source  level  1  (dB)  160 

Frequency  2 (hz)  300 

Source  level  2  (dB)  75 
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SPA  Data 

Latitude  50  00  N 

Longitude  020  00  W 

Semi-major  axis  (nm)  60.0 

Semi-minor  axis  (nm)  30.0 

Orientation  (deg  N)  90.0 

Containment  percent  86.0 

The  input  parameters  which  are  changed  to  generate  output  patterns  for  fleeing  da¬ 
tum  targets  is  presented  in  the  next  list. 

Setup 

Search  duration  (hrs)  8.0 

SPA  Data 

Semi-major  axis  (nm)  30.0 

Semi-minor  axis  (nm)  15.0 

The  two  input  parameters  that  are  changed  to  simulate  sonobuoys  with  detection 
ranges  which  van-  according  to  the  Lambda-Sigma  Jump  model  are  as  follows: 

Setup 

Ambient  acoustic  fluctuation  rate  (per  hour)  1.0 

Standard  deviation  of  ambient  acoustic  fluctuations  (dB)  8.0 
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APPENDIX  B.  SIMULATION  PROGRAMS 


-  This  appendix  presents  listings  of  examples  of  the  FORTRAN  computer  code  used 

in  the  simulation  of  the  sonobuoy  patterns  deployed  against  the  three  target  motion 
‘  models  of  SALT. 

The  two  programs  used  to  simulate  patrolling  target  motion  models  are  shown  here. 
The  first  program  shows  the  FORTRAN  code  which  tests  the  SALT  output  patterns 
generated  using  Definite  Range  Law  sonobuoys.  The  second  program  provides  a  model 
for  sonobuoys  which  behave  according  to  a  Lambda-Sigma  Jump  process.  The  pro¬ 
grams  used  to  simulate  the  transiting  target  and  the  fleeing  datum  target  are  similar  ex¬ 
cept  in  the  required  target  input  parameters. 

Following  the  patrolling  target  simulation  listings  is  the  listing  of  the  detection 
subroutine.  This  subroutine  is  used  by  all  six  of  the  simulation  programs  to  determine 
whether  or  not  the  sonobuoy  detection  criterion  is  met. 
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0 


1  PROGRAM  PATROL 

2  REAL  SIGMAX,  SIGMAY,  SCHTM,  LMBDA,  DETRNG 

3  INTEGER  SEED,  DET,  NUMDET,  I 

4  REAL  AVGSP ,  SPEED,  U,  COURSE,  XI,  Yl,  X2,  Y2,  XC,  YC 

5  REAL  T,  LEGTM,  TOTTIM,  DSTDET,  DT,  DETTIM 

6  CALL  EXCMSC'FILEDEF  01  DISK  PATPRM  DATA  A') 

7  CALL  EXCMSC'FILEDEF  02  DISK  PATBYS  DATA  A') 


8 

READ(01,*)  SIGMAX,  SIGMAY,  SCHTM,  LMBDA,  DETRNG,  SEED,  AVGSP 

9 

NUMDET  =  0 

10 

DO  40  I  =  I,  10000  , 

1 

11 

CALL  LNORMCSEED,  XI,  1,  1,  0) 

1 

12 

XI  =  XI  *  SIGMAX 

1 

13 

CALL  LNORMCSEED,  Yl,  1,  1,  0) 

1 

14 

Yl  =  Yl  *  SIGMAY 

1 

15 

TOTTIM  =  0.  0 

1 

16 

DETTIM  =  10000. 0 

1 

17 

10 

CONTINUE 

1 

18 

CALL  LRNDCSEED,  COURSE,  1,  1,  0) 

1 

19 

COURSE  =  COURSE  *  6. 283185 

1 

20 

CALL  LRNDCSEED,  U,  1,  1,  0) 

1 

21 

IF  (U  .LE.  1.  0/3.0)  THEN 

1 

1 

22 

SPEED  =  AVGSP*(24  +  SQRT(3024*U  +  16) )/56 

1 

1 

23 

ELSE 

1 

1 

24 

SPEED  =  AVGSP*(60  -  SQRT(1528  -  1512*U))/28 

1 

1 

25  : 

END  IF 

1 

26 

CALL  LEXPNCSEED,  T,  1,  1,  0) 

1 

27 

LEGTM  =  T/ LMBDA 

1 

28 

X2  =  XI  +  SPEED  *  LEGTM  *  SIN(COURSE) 

1 

29 

Y2  =  Yl  +  SPEED  *  LEGTM  *  COSCCOURSE) 

1 

30 

20 

CONTINUE 

1 

31 

READ(02,*)  XC,  YC 

1 

32 

IF  (XC  .GT.  999.0)  GO  TO  30 

1 

33 

CALL  DTCTC DETRNG,  XI,  Yl,  X2,  Y2,  XC,  YC,  DET,  DSTDET) 

1 

34 

IF  (DET  .EQ.  1)  THEN  ' 

1 

1 

35 

DT  =  TOTTIM  +  DSTDET/ SPEED 

1 

1 

36 

IF  (DT  .LT.  DETTIM)  THEN 

2 

1 

37 

DETTIM  =  DT 

2 

1 

38 

END  IF 

1 

1 

39 

END  IF 

1 

40 

GO  TO  20 

1 

41 

30 

CONTINUE 

1 

42 

REWIND  02 

1 

43 

XI  =  X2 

1 

44 

Yl  =  Y2 

1 

45 

TOTTIM  =  TOTTIM  +  LEGTM 

1 

46 

IF  (TOTTIM  .LT.  SCHTM  .AND.  DETTIM  . GT.  9999.0)  GOTO  10 

1 

47 

IF  (DETTIM  .LT.  SCHTM)  THEN 

1 

1 

48 

NUMDET  =  NUMDET  +  1 

1 

1 

49 

WRITE( 11 ,*)  DETTIM 

1 

1 

50 

END  IF 

1 

51 

40 

CONTINUE 

52 

PRINT*,  'NUMBER  OF  DETECTIONS  =’ ,  NUMDET 

53 

STOP 

54 

END 

3R 


1 

PROGRAM  LSPAT 

2 

REAL  SIGMAX,  SIGMAY,  SCHTM,  LMBDA,  DETRNG ,  FOM 

3 

REAL  ACFL,  AFTIM,  R,  Rl,  TL,  TL1 

4 

INTEGER  SEED,  DET,  NUMDET,  I 

5 

REAL  AVGSP,  SPEED,  U,  COURSE,  XI,  Yl,  X2,  Y2, 

XC,  YC 

6 

REAL  T,  LEGTM ,  TOTTIM,  DSTDET,  DT,  DETTIM 

7 

CALL  EXCMSC'FILEDEF  01  DISK  PATPRM  DATA  A1) 
CALL  EXCMSC'FILEDEF  02  DISK  PATBYS  DATA  A') 

8 

9 

CALL  EXCMSC'FILEDEF  03  DISK  TLCRV  DATA  A') 

10 

READCOl,*)  SIGMAX,  SIGMAY,  SCHTM,  LMBDA,  FOM, 

SEED,  AVGSP 

11 

NUMDET  =  0 

12 

DO  110  1=1,  10000 

1 

13 

CALL  LNORMCSEED,  XI,  1,  1,  0) 

1 

14 

XI  =  XI  *  SIGMAX 

1 

15 

CALL  LNORMCSEED,  Yl,  1,  1,  0) 

1 

16 

Yl  =  Yl  *  SIGMAY 

1 

17 

TOTTIM  =  0.  0 

1 

18 

DETTIM  =  10000. 0 

1 

19 

CALL  LRND(SEED,  COURSE,  1,  1,  0) 

1 

20 

COURSE  =  COURSE  *  6. 283185 

1 

21 

CALL  LRND(SEED,  U,  1,  1,  0) 

1 

22 

IF  (U  .LE.  1.0/3.  0)  THEN 

1 

1 

23 

SPEED  =  AVGSP*( 24  +  SQRT(3024*U  +  16))/56 

1 

1 

24 

ELSE 

1 

1 

25 

SPEED  =  AVGSP* ( 60  -  SQRTC 1528  -  1512*U))/28 

1 

1 

26 

END  IF 

1 

27 

CALL  LEXPNCSEED,  T,  1,  1,  0) 

1 

28 

LEGTM  =  T/LMBDA 

1 

29 

CALL  LEXPNCSEED,  AFTIM,  1,  1,  0) 

1 

30 

10 

CALL  LNORMCSEED,  ACFL,  1,  1,  0) 

1 

31 

ACFL  =  ACFL  *  8.  0 

1 

32 

FOM  =  FOM  +  ACFL 

1 

33 

IF  (FOM  .LE.  66.0  .OR.  FOM  .GE.  100.66)  GO 

TO  10 

1 

34 

20 

CONTINUE 

1 

35 

READC03,*)  R,  TL 

1 

36 

IF  (FOM  .GE.  TL)  THEN 

1 

1 

37 

Rl  =  R 

1 

1 

38 

TL1  =  TL 

1 

1 

39 

ELSE 

1 

1 

40 

DETRNG  =  Rl  +  (FOM  -  TL1)/(TL  -  TL1) 

1 

1 

41 

GO  TO  30 

1 

1 

42 

END  IF 

1 

43 

GO  TO  20 

1 

44 

30 

CONTINUE 

1 

45 

REWIND  03 

1 

46 

35 

CONTINUE 

1 

47 

IF  (AFTIM  .LT.  LEGTM)  THEN 

1 

1 

48 

X2  =  XI  +  SPEED  *  AFTIM  *  SIN(COURSE) 

1 

1 

49 

Y2  =  Yl  +  SPEED  *  AFTIM  *  COS(COURSE) 

1 

1 

50 

40 

CONTINUE 

1 

1 

51 

READC02,*)  XC,  YC 

1 

1 

52 

IF  CXC  .GT.  999. 0)  GO  TO  50 

1 

1 

53 

CALL  DTCT( DETRNG,  XI,  Yl,  X2,  Y2,  XC, 

YC,  DET,  DSTDET) 

1 

1 

54 

IF  (DET  .EQ.  1)  THEN 

2 

1 

55 

DT  =  TOTTIM  +  DSTDET/SPEED 

2 

1 

56 

IF  (DT  .LT.  DETTIM)  THEN 

39 


3 

1 

57 

DETTIM  =  DT 

3 

1 

58 

END  IF 

2 

1 

59 

END  IF 

1 

1 

60 

GO  TO  40 

1 

1 

61 

50 

CONTINUE 

1 

1 

62 

REWIND  02 

1 

1 

63 

XI  =  X2 

1 

1 

64 

Y1  =  Y2 

1 

1 

65 

LEGTM  =  LEGTM  -  AFTIM 

1 

1 

66 

TOTTIM  =  TOTTIM  +  AFTIM 

) 

1 

1 

67 

CALL  LEXPN( SEED ,  AFTIM,  1,  1,  0) 

1 

1 

68 

60 

CALL  LNORM(SEED,  ACFL,  1,  1,  0) 

1 

1 

69 

ACFL  =  ACFL  *  8.  0 

1 

1 

70 

FOM  =  FOM  +  ACFL 

1 

1 

71 

IF  (FOM  .  LE.  66. 0  . OR.  FOM  .  GE.  100. 66)  GO  TO  60 

t 

1 

1 

72 

70 

CONTINUE 

1 

1 

73 

READ(03,*)  R,  TL 

1 

1 

74 

IF  (FOM  .GE.  TL)  THEN 

2 

1 

75 

R1  =  R 

2 

1 

76 

TL1  =  TL 

2 

1 

77 

ELSE 

2 

1 

78 

DETRNG  =  R1  +  (FOM  -  TL1)/(TL  -  TL1) 

2 

1 

79 

CO  TO  80 

2 

1 

80 

END  IF 

1 

1 

81 

GO  TO  70 

1 

1 

82 

80 

CONTINUE 

1 

1 

83 

REWIND  03 

1 

1 

84 

ELSE 

1 

1 

85 

X2  =  XI  +  SPEED  *  LEGTM  *  SIN(COURSE) 

1 

1 

86 

Y2  =  Yl  +  SPEED  *  LEGTM  *  COS( COURSE) 

V 

1 

1 

87 

90 

CONTINUE 

1 

1 

88 

READ(02,*)  XC,  YC 

1 

1 

89 

IF  (XC  .GT.  999.0)  GO  TO  100 

dstiSet 

1 

1 

90 

CALL  DTCT( DETRNG,  XI,  Yl,  X2,  Y2,  XC,  YC,  DET, 

1 

1 

91 

IF  (DET  .EQ.  1)  THEN 

2 

1 

92 

DT  =  TOTTIM  +  DSTDET/SPEED 

2 

1 

93 

IF  (DT  .LT.  DETTIM)  THEN 

3 

1 

94 

DETTIM  =  DT 

3 

1 

95 

END  IF 

2 

1 

96 

END  IF 

1 

1 

97 

GO  TO  90 

1 

1 

98 

100 

CONTINUE 

1 

1 

99 

REWIND  02 

1 

1 

100 

XI  =  X2 

1 

1 

101 

Yl  =  Y2 

1 

1 

102 

AFTIM  =  AFTIM  -  LEGTM 

1 

1 

103 

TOTTIM  =  TOTTIM  +  LEGTM 

1 

1 

104 

CALL  LRND(SEED,  COURSE,  1,  1,  0) 

1 

1 

105 

COURSE  =  COURSE  *  6. 283185 

1 

1 

106 

CALL  LRND( SEED ,  U,  1,  1,  0) 

1 

1 

107 

IF  (U  .LE.  1.0/3.  0)  THEN 

2 

1 

108 

SPEED  =  AVGSP*(24  +  SQRT(3024*U  +  16))/56 

* 

2 

1 

109 

ELSE 

2 

1 

110 

SPEED  =  AVGSP*(60  -  SQRT(1528  -  1512*U))/28 

2 

1 

111 

END  IF 

/ 

1 

1 

112 

CALL  LEXPN(SEED,  T,  1,  1,  0) 
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113  LEGTM  =  T/LMBDA 

114  END  IF 

115  IF  (TOTTIM  .LT.  SCHTM  .AND.  DETTIM  .  GT.  9999.0)  GO  TO  35 

116  IF  (DETTIM  .  LT.  SCHTM)  THEN 

117  NUMDET  =  NUMDET  +  1 

118  WRITE( 11  ,*)  DETTIM 

119  END  IF 

120  110  CONTINUE 

121  PRINT* ,  'NUMBER  OF  DETECTIONS  =' ,  NUMDET 

122  STOP 
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11 
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16 
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18 

19 

20 
21 
22 

23 

24 

25 
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27 

28 

29 

30 

31 

32 

33 

34 
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36 

37 


X2)  THEN 


,EQ.  Y2)  THEN 


SUBROUTINE  DTCT(R,  XI,  Yl,  X2,  Y2,  XC,  YC,  DET,  DSTDET) 

REAL  R,  RSQ,  XI,  Yl,  X2,  Y2,  XC,  YC,  XP,  YP,  M,  DSQ 
REAL  D12SQ,  D1PSQ,  D2PSQ,  Z,  DSTDET 
INTEGER  DET 
RSQ  =  R  *  R 
IF  (XI  .EQ. 

XP  =  XI 
YP  =  YC 
ELSE  IF  (Yl 
XP  =  XC 
YP  =  Yl 
ELSE 

M  =  (Y2  -  Y1)/(X2  -  XI) 

XP  =  (XC/M  +  YC  +  M*X1  -  Y1)/(M  +  1/M) 

YP  =  M*XP  -  M*X1  +  Yl 
END  IF 

DSQ  =  (XC  -  XP)*(XC-XP)  +  (YC  -  YP)*(YC  -  YP) 

DET  =  0 

IF  (DSQ  .  LT.  RSQ)  THEN 
D12SQ  =  (XI  -  X2)*(X1  - 
D1PSQ  =  (XI  -  XP)*(X1  - 
D2PSQ  =  (XP  -  X2)*(XP  - 
IF  ((XC  -  X1)*(XC  -  XI)  +  (YC  -  Y1)*(YC  -  Yl)  .  LT.  RSQ)  THE: 
DET  =  1 
DSTDET  =  0.  0 

ELSE  IF  ((D1PSQ  .  LT.  D12SQ) 

DET  =  1 

Z  =  SQRT(RSQ  -  DSQ) 

DSTDET  =  SQRT(DIPSQ) 

ELSE  IF  ((X2  -  XC)*(X2 
l-  THEN 

DET  =  1 

Z  =  SQRT(RSQ  -  DSQ) 

DSTDET  =  SQRT(DIPSQ) 

END  IF 
END  IF 
RETURN 
END 


X2) 

XP) 

X2) 


(Yl  -  Y2)*(Y1  -  Y2) 
(Yl  -  YP)*( Yl  -  YP) 
(YP  -  Y2)*(YP  -  Y2) 


-  Z 

XC) 


-  Z 


AND.  (D2PSQ  .LT.  D12SQ))  THEN 


+  (Y2  -  YC)*( Y2  -  YC)  . LT.  I^sd 
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